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Description 

The present invention relates to systems for pro* 
viding, by nuclear magnetic resonance, images of 
distributions of a chosen quantity in a selected 
region of a body. The examination may be of 
many different kinds of bodies but a particularly 
beneficial application, to which much attention 
has been given, is medical examination of pa- 
tients. 

For NMR examination, suitable combinations 
of magnetic fields are applied to the body to be 
examined by appropriate magnet (coil) systems. 
Excited material is then detected by currents in- 
duced in one or more detector coil systems which 
are analysed to provide the required distribution. 

Various schemes have been proposed for this 
purpose, two of which are relevant to this inven- 
tion being that proposed by Edelstein and Hutchi- 
son et al In "Physics in Medicine & Biology" 25, 
751 - 756 (1980) and that invented by Ernst and 
described in U.S. Patent No. 4 070 61 1 . 

The Imaging sequence used for the spin warp 
NMR imaging proposed by Edelstein et al shown 
in Figure 1 produces spin-echo images with 
transaxial orientation. The sequence begins when 
a sinc-modulated 90° nutation pulse (A) occurs 
simultaneously with a Gz gradient pulse to pro- 
duce the spatially selective excitation of nuclear 
spins. Three gradient pulses (B) are applied; Gy 
phase encodes in-plane resolution, G z phase en* 
codes Image section, and G x defocuses the free 
induction decay (FID) to temporarily centre the 
later occurring spin-echo signal. A sinc-modu- 
lated 180° RF pulse (C) occurs simultaneously 
with a Gz gradient pulse to produce echo refocus- 
ing. A read-out gradient G x frequency encodes 
(D) the spin-echo signal for resolution along an 
inplane axis. The 2 D images were made by itera- 
tion of the pulse sequence each time with a differ- 
ent magnitude of Gy at (B). 

Each of the above described cycles of the 
pulse sequence can be repeated several times 
(N) for signal averaging and to cancel FID arti- 
facts which are produced by FID signal active 
during the spin-echo but generated by the 90° 
irradiation of the planar region outside the inner 
volume. To produce sagittal images, ail the gra- 
dient pulses Gz and Gx are interchanged. 

It should be appreciated that the above de- 
scribed techniques are disadvantageous in vari- 
ous aspects. For example, In the spin warp NMR 
imaging or Fourier zeugmatography known from 
U.S. Patent No. 4 070 61 1 two phase encoding 
magnetic gradients must be varied linearly in the 
respective directions so as to obtain the Image 
Information depending upon the linear variation 
between the intensity of the magnetic fields and 
the corresponding coordinates. However, in prac- 
tice, the linearity depends on various instrumental 
parameters, such as quality of the coil, the power 
supply for the coil, etc. Therefore, the intensity of 
the gradient fields does not vary in a linear man- 
ner around the original which is the point in which 
the gradient contributes to the external field. This 



unavoidable non linearity of the gradient causes 
artifacts in the reconstructed NMR images. 

It is an object of this invention to provide an 
improved method for producing nuclear magnetic 
5 resonance image data. 

More specifically, an object of the present in- 
vention is to provide a method which eliminates 
those artifacts due to the above described non-li- 
nearity in the gradient, and includes a correction 

10 factor to correct variation of the intensity of the 
magnetic gradient along one axis prior to collec- 
tion of the signals required for forming an image 
of magnetisation distribution using the Fourier 
zeugmatography technique. 

15 According to the invention, prior to collecting 

an echo resonance signal produced by the 180° 
pulse, when the so-called spinecho technique is 
used, which is subject to 2 D Fourier transforma- 
tion to obtain an image of the magnetisation dis- 

20 tribution in the selected slice, a preliminary modu- 
lated phase encoding gradient field is applied, so 
that several preliminary and reference echo sig- 
nals are obtained by sampling operations being 
performed in the presence of the preliminary gra- 

25 dient magnetic fields. Then, the intensity of the 
gradient with reference to the obtained signal is 
used to start a phase encoding gradient pulse 
sequence therewith for collecting the echo sig- 
nals actually used for forming the desired im- 

30 aging. 

One embodiment of the invention will now be 
described by way of example with reference to 
the accompanying drawings, in which: 

35 Figure 1 shows the pulse sequence of a spin 
warp NMR imaging method; 
Figure 2 diagrammatical ly illustrates the ar- 
rangement of a target within the mag- 
netic field generating coils; 

40 Figure 3 shows schematically an apparatus 
which can be used to implement the 
invention; 

Figure 4 shows the pulse sequence and spin- 
echo signals provided by the method 
45 of this invention; and 

Figure 5 shows the relationship between the 
practical amplitude of the phase-en- 
coding gradient, and range of ampli- 
tude of the gradient used in acoord- 
50 ance with this invention. 

Referring to Figure 2, a target or object to be 
investigated is placed in a static magnetic field Bo 
which extends in the 2-directIon. RF transmitter 

55 coils 4 and 6 and RF receiver coils 8 and 10 are 
shown, for example, to consist of two separate 
coil sets that are electrically orthogonal to one 
another and which generate RF magnetic field 
pulses having the required modulation to excite 

60 resonance in a planar slice 12 of the object 
defined by the combination of the application of 
the magnetic field gradient Gz along the Z-direc- 
tion produced by Z-axis gradient coils 16a and 
16b and application of the corresponding RF 

65 magnetic field. RF magnetic field pulses are di- 
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rected orthogonal to the Bo field. The Z-axis gra- 
dient coils 16a and 16b are in the form of a Helm- 
holtz pair wound in opposition, and are disposed 
substantially jn planes perpendicular to the Z 
axis. In addition to G z gradient coils 16a and 16b, 
another two gradient coil sets are necessary to 
provide gradient G x . and Gy, respectively for en- 
coding the phase of the spin echo signals for 
image reconstruction to spatially resolve informa- 
tion along each axis. For the sake of clarity, the 
known arrangement of Gx and Gy gradient coil 
sets is only shown schematically as 16c f 16d. 
16e, and 161 respectively in Figure 3, but Is 
shown.in more detail for example, in U.S. Patent 
4 355 282. 

Figure 3 shows an apparatus for performing 
NMR detection in accordance with one embodi- 
ment of the invention. The apparatus comprises 
coils 16c. 16d, 16e ( and 16f gradient field driver 
18 for the sexcitation of the coils 16a. 16b. 16c. 
16d f 16e and 16f, respectively, under the control 
of sequence controller 20, and a high-frequency 
signal generator 22 for the excitation of the coils 
4, 6. The apparatus also comprises a high-fre- 
quency signal detector 24, processing means 26 
Including a sample circuit, an analogue-to-digital 
converter, a memory, and a Fourier transformer 
for performing a Fourier transform atbn. The static 
magnetic field B 0 which has already been de- 
scribed with reference to Figure 2 is generated by 
a static magnetic field generator 30 including ap- 
propriate coils and electrically connected driver. 
Though only the G x . Gy gradient coils 16c. 16d, 
16e, and 16f are Illustrated schematically In Fig- 
ure 3, the G x coils 16c and 16d. and the Gy coils 
16e and 16f are two sets of coils axially dis- 
placed, each set comprising two pairs of saddle 
coils the pair (16c, 16d) being at 90° to the pair 
(16e, 161). These coils are themselves inserted 
into the central aperture in B 0 coils (not shown in 
Figures). RF coil sets 4, 6, 8 and 10 are two of 
saddle shaped coils which are driven In parallel, 
to provide the rotating RF field by the highfre- 
quency signal generator 22. and also to detect 
the echo resonance signals, which are of approxi- 
mately the same frequency, by the highfrequency 
signal detector 24. 

The apparatus performs a method for deter- 
mining the nuclear magnetisation distribution in a 
body 2 as will be described hereinafter. The 
method comprises several steps. 

The pulse sequence used to form an image 
from single spin echos following separate excita- 
tions Is shown In Figure 4. 



Interval 1 

A selective 90° RF pulse is applied simultaneous- 
ly with a magnetic field gradient G 2 . This selec- 
tively excites nuclear spins In and close to the 
plane Z - Zo. The value of Z> can be altered by a 
change in the frequency of the 90° pulse. 



4 

Interval 2 

A gradient of the magnetic field having a positive 
value G x is applied to dephase the nuclear spins 
5 along the X-direction. Simultaneously a smaller 
signal gradient Gy(p) of the magnetic field is ap- 
plied to dephase the spins along Y as a prelimi- 
nary prior to collection of the resonance signals 
used for forming an image. 

10 

Interval 3 

A 180° RF pulse Is applied simultaneously with a 
15 longer Gz gradient pulse to produce echo refo- 
cusing. During a part of the period after the 180° 
pulse, a magnetic field gradient Gy*(p) having 
corresponding values of time integral of the gra- 
dient Gy(p) applied before the 1 80° pulse is ap- 
20 plied. 



Interval 4 

25 The echo signal is measured and sampled in the 
presence of a readout gradient G x which is, ap- 
plied to cause the spins to rephase and produce 
a kind of gradient-induced spin echo in the middle 
of the signal acquisition interval. The Fourier 

30 transform of this spin echo signal is therefore a 
one-dimensional projection of the spin density 
within the slice on the the X axis. 

The entire set of steps described above is 
successively repeated a number of times and In 

35 each repetition a different amplitude of gradient 
Gy'(p) is utilised to cause the signal to reach a 
peak value, of which amplitude represents the 
designated (nominal) phase encoding of the sig- 
nal to spatially resolve information along the Y 

40 axis. 



Interval 5 

45 System recovery time until the occurrence of In- 
terval 1 of the next sequence. This should be long 
compared with the spin-lattice relaxation time T 
and may be of the order of a second in a whole- 
body imaging apparatus. 

50 The same procedures are performed to the 

reversed magnetic gradients Gy(p) and Gy'(p). 

The resultant gradient Gy(p) max serves to 
give a desired amount of twist or "warp" to each 
vertical column of spins (Y axis vertical) defined 

55 by the presence of the gradient G x . In fact it maxi- 
mises the response to a particular vertical spatial 
frequency in the column equal to the spatial fre- 
quency of the "warp*. 

After the above described preliminary step, in 

60 Figures 4 and 5 the phase-encoding gradient Gy 
is always applied for the same period of time but 
changes in strength for different pulse sequen- 
ces. In each successive sequence Gy has the 
same shape and length but Hs amplitude changes 

65 by equal steps from an imaginary zero, where the 
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effect of two phase encoding gradients G y (p) and 
Gy'(p) Is cancelled in phase of the echo signal, to 
a maximum value. 

In accordance with the proposed method, a 
range of vertical frequencies from zero up to a 
maximum is a linear function of G y linear phase- 
encoding gradient strengths starting with the 
value of G/(p) max by avoiding the non-linear 
scale of phase-encoding gradient. 

The two-dimensional (2D) Fourier transform 
image is produced by iteration of the pulse se- 
quence each time with a different magnitude of 

In order to achieve this object of applying G y 
linear phaseencoding gradient strengths starting 
with the point where the phase of the encoded 
spin-echo signal is of imaginary zero by the com- 
bination of application of preliminary gradients 
G y (p) before and after 180° RF pulse, the spin- 
echo signals are detected by the high-frequency 
detector 24 In order of changing G y '(p) pulse 
strength generated by Gy coils 16c and 16f under 
control of sequence controller 20. Then, In pro- 
cessing means (26 in Figure 3) the time Integral 
value of G/(p) gradient is determined, which 
maximises the signal that Is eventually obtained. 
Information on this value of the phase-encoding 
gradient is provided to the sequence controller 20 
to follow the primary image data acquisition steps 
with application of different Gy linear phaseencod- 
ing gradients from the offset magnetic field that 
has been previously obtained. 



Claim 

Method of deriving image Information from an ob- 
ject using nuclear magnetic resonance signals 
comprising subjecting said object to a static mag- 
netic field along an axis and carrying out a set of 
sequential preliminary steps which comprises: 

- selectively exciting nuclear spins in a selected 
plane In the presence of a first gradient ma- 
gnetic field (G*) applied in a first direction; 

- applying a second gradient magnetic field for 
(Gx) in a second direction encoding the phase 
of the selected spins to spatially resolve infor- 
mation along one axis which Is orthogonal to 
said first direction; 

- applying a preliminary third gradient magnetic 
field (Gy) In a third direction orthogonal to said 
first and second directions for preliminary 
phase encoding along said third direction; 

- applying a 180° RF pulse associated with a 
gradient applied in the first direction to 
produce echo refocussing; 

- applying an additional preliminary third gra- 
dient (GV(p)) after said 180° RF pulse having 
a designated strength which Is substantially 
identical to that of said preceding preliminary 
third gradient (Gy); 

- removing the additional third gradient (GV(p)) 
and reading out the resultant spin echo signal 
from the object In the presence of a read-out 



6 

gradient applied in the second direction; and 
then 

- successively repeating the above set of se- 
quential steps with a different value of Integral 

5 of said additional preliminary third gradient 

(G'y(p)) applied after the 180° RF pulse in 
order to calculate on the basis of a processing 
of resultant spin echo signals an intensity of 
the phaseencoding gradient in the third dlreo 

10 tlon preserving Its linearity and then successi- 

vely repeating the above set of sequential 
steps with intensities of the preliminary third 
gradient replaced by different linear phase-en- 
coding gradient intensities starting with said 

15 obtained Intensity to obtain a group of 

sampled signals wherefrom, after Fourier 
transformation thereof, an image of the distri- 
bution of the spin magnetisation is determin- 
ed. 

20 

, Patentanspruch 

Verfahren zur Gewinnung von Bildinformationen 
25 elnes Objektes aus magnetischen Kemresonanz- 
signalen, Indem das Objekt einem statischen 
Magnetfeld entlang elner Achse ausgesetzt wlrd 
und eine Reihe von aufelnanderfolgenden vor- 
bereitenden Schritten durchgefOhrt werden, die 
30 folgende Schritte umfassen: 

- Selective Erregung von Kernspins in einer 
ausgewaMen Ebene bei Vorhandensein 
eines ersten, in einer ersten Richtung ange- 

35 legten Magnetfeldgradienten (G z ); 

- Anlegen eines zweiten Magnetfeldgradienten 
(Gz) In einer zweiten Richtung zum Kodieren 
der Phase der ausgewahlten Spins, urn die 
Information raumlich aufzulfisen entlang einer 

40 Achse, die orthogonal zur ersten Richtung 

verlauft; 

- Anlegen eines vorbereitenden dritten Magnet- 
feldgradienten (Gy) in einer dritten Richtung, 
die orthogonal zur ersten und zweiten Rich- 

45 tung verlauft, zur vorbereitenden Phasenko- 

dierung entlang der dritten Richtung; 

- Anlegen eines 180° Hochfrequenzimpulses, 
der einem In der ersten Richtung angelegten 
Gradienten zugeordnet 1st, urn eine Echo-Re- 

50 fokussierung zu erzeugen; 

- Anlegen eines zusatzllchen vorbereitenden 
dritten Gradienten (G y(p)) nach Anlegen des 
180° Hochfrequenzimpulses, mit einer be- 
stlmmten Starke, die Im wesentlichen gleich 

55 1st derjenigen des vorhergehenden, vorberei- 

tenden dritten Gradienten (Gy); 

- Beseitigung des zusatzlichen dritten Gradien- 
ten (Gy(p)) und Auslesen des resuhlerenden 
Spin-Echo-SignaJes aus dem Objekt bei Vor- 

60 handenseln eines in der zweiten Richtung an- 

gelegten Auslese-Gradienten; und 

- Wiederholung der Reihe nach des oberen 
Satzes aufeinanderfolgender Schritte mit 
einem IntegraJwert des zusatzlichen nach 

65 dem 180° Hochfrequenzimpufs angelegten 
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Gradients n (G' y (p)), urn auf der Grundlage 
einer resultierenden Spin-Echo-Signalverar- 
beitung erne Intensitat des phasenkodieren- 
den Gradienten In der dritten, seine Llnearitat 
erhaltenden Rlchtung zu berechnen, und 
dann aufeinanderfolgende Wiederholung des 
oberen Satzes aufeinanderfolgender Schritte, 
wobel die Intensitaten des vorbereitenden 
dritten Gradienten ersetzt sind durch unter- 
schiedliche lineare phasenkodierende Gra- 
dienten-lntensilaten, indem mit der erhaltenen 
Intensitat begonnen wird, um eine Gruppe 
von abgetasteten Signalen zu erhalten, aus 
der nach Fourier-Transformation ein Vertei- 
lungsbild der Spln-Magnetislerung ermhtell 
wird. 



8 

signaux de spins resultant, une intensite du 
gradient de codage de phase dans la troisie- 
me direction en preservant sa linearite, puis a 
repeter successivement Fensemble cl-dessus 

5 d'etapes sequentielles avec des intensites du 

troisieme gradient preliminaire remplacees 
par differentes intensites de gradients de 
codage de phase lineaires en commencant 
par ladite intensite obtenue, afin d'obtenir un 

10 groupe de signaux echantillonnes, a partir 

desquels, et apres leur transformation par une 
transformee de Fourier, on determine la re- 
partition de la magnetisation du spin. 

15 



Revendicatlon 

20 

Un proc€de destine a obtenir des informations 
d'images a partir d'un objet en utilisant des sig- 
naux de resonance magnetique nucleaire, en 
soumettant ledit objet a un champ magnetique 
statique le tang d'un axe et en effectuant un en- 25 
semble d'etapes preliminaires sequentielles qui 
consistent : 

- a excher de f agon selective des spins nucleai- 

res dans un plan cnoisi en presence d'un 30 
premier champ magnetique de gradient Gi 
applique dans une premiere direction; 

- a appliquer un second champ magnetique de 
gradient G x dans une seconde direction pour 

coder ia phase des spins choisls afin de tra- 35 
duire des informations dans Tespace le long 
d'un axe qui est orthogonal a ladite premiere 
direction; 

- a appliquer un troisieme champ magnetique 

de gradient preliminaire Gy dans une troisie- 40 
me direction orthogonale auxdttes premiere et 
seconde directions pour un codage en phase 
preliminaire le long de ladite troisieme direc- 
tion; 

- a appliquer une impulsion a haute frequence 45 
a 180° associee a un gradient applique a la 
premiere direction pour produire une refocali- 
sation cfecho; 

- appliquer un troisieme gradient preliminaire 
supplementaire (G' y (P)) apres ladrte impulsion 50 
a haute frequence a 180°, ayant une Intensite 
definie qui est sensiblement identique a cede 

du-dit troisieme gradient precedent preliminai- 
re (Gy); 

- a supprimer le troisieme gradient preliminaire 55 
(G'y(P)) et a extraire le signal d'echo de spin 
resultant a partir de Tobjet en presence d'un 
gradient extra'rt applique dans une seconde 
direction; puts. 

- a repeter successivement fensembie ct-des- 60 
sus d'etapes sequentielles avec une valeur 
dtfferente d'integrale du-dit troisieme gradient 
preliminaire supplementaire (G'y(p)) applique 

apres Hmpulsion a haute frequence a 180° en 

vue de calculer sur la base d'un trartement de 65 
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Fig. I. 
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